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Abstract 

In this work it is investigated the evolution of a Universe where a scalar field, non-minimally 
coupled to space-time curvature, plays the role of quintessence and drives the Universe to 
a present accelerated expansion. A non-relativistic dark matter constituent that interacts 
directly with dark energy is also considered, where the dark matter particle mass is assumed 
to be proportional to the value of the scalar field. Two models for dark matter pressure are 
considered: the usual one, pressureless, and another that comes from a thermodynamic theory 
and relates the pressure with the coupling between the scalar field and the curvature scalar. 
Although the model has a strong dependence on the initial conditions, it is shown that the 
mixture consisted of dark components plus baryonic matter and radiation can reproduce the 
expected red-shift behavior of the deceleration parameter, density parameters and luminosity 
distance. 

1 Introduction 

The measurements of the rotation curves of spiral galaxies as well as other astronomical experi- 
ments suggest that the luminous matter represents only a small amount of the massive particles of 
the Universe, and that the more significant amount is related to dark matter. 

Recently the astronomical observations with super-novae of type la suggested that our Universe 
is presently submitted to an accelerated expansion [21 El; the nature of the responsible entity, called 
dark energy, still remains unknown. The simplest explanation for the acceleration is a cosmological 
constant (see 0]), which fits the present data very well but has some important unsolved problems. 
Another possibility is to introduce a scalar field <j>(t), which has been extensively studied by the 
scientific community. By considering a barotropic equation of state for the scalar field = upj,, 
we have a constant value of lu = —1 for the cosmological constant model and a variable lu > — 1 for 
a minimally coupled scalar field model. However, the measured data from Hubble Space Telescope 
[S] states the restriction for lu = — 1.02+q'jJ, whereas some recent observations jS] has drawn 
the attention to the viability of models where lu < — 1, which would invalidate both cosmological 
constant and minimally coupled scalar field descriptions. Some phantom field models were proposed 
in order to contemplate this hypothesis [7j, but they face also some strong difficulties 

In this work we consider a scalar field non-minimally coupled to space-time curvature, which 
was investigated in and widely studied recently, among others, in the works ^3 111! 1121 1131 ITi] . 
In a more recent paper |15| the suitability of such models was studied and it was shown that the 
Einstein field equations are stable under some conditions for the constant that couples the scalar 
field with the curvature scalar. 

Aside several models for the dark sector interaction 1171 1181 ITU) - that could replace the 
cosmological constant model in the case of future experimental tests - we consider here a direct 
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coupling between dark matter and dark energy. In this model, studied in |2(J| and more recently in 
|21j . the dark matter particle mass is proportional to the value of the scalar field that represents 
the dark energy. However, in this model it is necessary to have some precautions, since by assuming 
the dark matter particle mass to change with time, we could make the dark matter energy density 
to become physically inconsistent at early stages. We emphasize that in this work the full set 
of differential equations are solved numerically, unlike the asymptotic approximation of the work 
|21j . The baryonic matter and radiation are described here as usual non- interacting components, 
being the first non-relativistic and pressureless. However, for the dark matter two models arc 
considered: the first one, by assuming the dark matter as pressureless, and the second one, by 
using a thermodynamic theory in order to relate the effects of the non-minimally coupling to the 
dark matter pressure. All components will be described by a set of field equations, and the resulting 
observables, - i.e., the density parameters, the decelerating parameter and the luminosity distance 
- which are obtained as solutions of the field equations will be compared to the available data set in 
order to drawn the conclusions about the viability of this model. We show that physically acceptable 
solutions are obtained, if we chose initial conditions restricting the rolling nature of the scalar field 
at early stages. Furthermore, there exist some freedom parameters that will be important to fit 
the data from incoming experiments. The work is organized as follows: in section 2 the basic field 
equations for each component, i.e. dark energy, dark matter, baryons and radiation, are derived. 
The solutions for the pressureless dark matter are found in section 3 for given initial conditions, 
coupling and potential constants, whereas in section 4 we give the corresponding solutions for the 
dark matter with non-vanishing pressure. We close the work with section 5 where we address to 
some remarks and summarize the results of previous sections. Units have been chosen so that 
8ttG = c = H = k = 1, whereas the metric tensor has signature (+,—,—,—). 



2 Field equations 

We shall consider that the Universe is modeled as a mixture of a scalar field which plays the role of 
dark energy, a dark matter field and non-interacting baryons and radiation. Here we adopt a more 
general scalar field which is non-minimally coupled to curvature [5] , and a dark matter field whose 
particle mass is proportional to the value of the scalar field at each point The action for this 
model is written as: 



{ \{i - ^ 2 )R + \d^d^ -v{4>) + \d^x - \c 2 < 

where <fi and x are the dark energy and dark matter scalar fields, respectively, V (</>) is the dark 
energy potential density, £ is the coupling constant between dark energy and the curvature scalar R 
and C,4> is the dark matter particle mass, with Q being a constant. Furthermore, C r and Cb denote 
the lagrangian densities of the radiation and of the baryons, respectively. 

By taking the variation of the action Q with respect to the metric tensor g^ Ul it follows Einstein 
field equations 

1 T „ 

R ^-- 2 R 9^ = - {1 _ V m ( 2 ) 

In the above equation the total energy-momentum tensor 7j,„ of the sources of the gravitational 
field is a sum of the energy- momentum tensors of the baryons T^„, radiation T* v , dark matter T^J 1 
and scalar field T*, i.e., = T^ v + T^ v + + T$. The expressions for the energy- momentum 
tensors read 

-,b 2 5y/=g£ b 2 5^f=g~L r 

T*? = d^X - [\doXd°X - l( 2 > 
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We call attention to the fact that due to the inclusion of the coupling between the scalar field 
and the scalar curvature, the covariant divergence of the total energy-momentum tensor does not 
vanish. Indeed, one can obtain from J5J by using Bianchi identities 

V V T^ = 2& (sr - \ng^ d v <t>, (6) 

indicating that the total energy-momentum tensor is not a conserved quantity. However, we shall 
assume that the covariant divergence of the energy-momentum tensors of the baryons and radiation 
vanish, i.e., V ' JT^ V = and V ' JT^ iV = 0, since they are considered as non- interacting fields from 
the decoupling age up to the present time. 

The evolution equations for the scalar and dark matter fields are obtained from Euler-Lagrangc 
equations, yielding 

v CT v ff + ^ + eR0 + CW = o, v CT v ff x + CV 2 x = o. (7) 

Henceforth, we shall restrict ourselves to the study of the consequences of the above model on 
a Universe which is isotropic, homogeneous and spatially flat, described by the Robertson- Walker 
metric: 

ds 2 = dt 2 - a{t) 2 (dx 2 + dy 2 + dz 2 ) , (8) 

where a (t) denotes the cosmic scale factor. 

For an isotropic and homogeneous Universe the most general representation for the total energy- 
momentum tensor in a comoving frame is given by {T^ v ) = diag(p, — p, — p, — p), where the total 
energy density p and pressure p of the mixture are given in terms of the respective quantities for 
its constituents, namely, p = Pb + Pr + Pdm + P<t> an( l P — Pb + Pr + Pdm + P<f>- From now on, we 
shall assume that the baryons are non-relativistic particles so that pb = and that the barotropic 
equation of state for the radiation field p r = p r /3 holds. 

Now for the Robertson- Walker metric, Einstein field equations (0) lead to the following modified 
forms of the Friedmann and acceleration equations 

H 2 = 9 - = P + 3P (Q) 

3(1-^ 2 )' a 6(1-^ 2 )' {1 

respectively, due to the non-minimally coupling of the scalar field to the gravitational field. Above 
H = d(t)/a(t) is the Hubble parameter and the dot denotes a differentiation with respect to time. 

The hypothesis of homogeneity implies that the fields (j> and \ must be only functions of time 
so that the evolution equations for the dark energy Qi and dark matter 02 fields reduce to 

4> + 3Hj> + ^ + (p - 3p) + i ( Pdm - Pdm ) - 0, (10) 

X + SHx + CVx = o, (11) 

respectively. 

The energy densities and the pressures of the dark energy and dark matter fields can now be 
determined from an d ©, yielding 

W = ^IHW, W -i0 2 -V-2£(# + 2 + 2w), (12) 

Pdm = \x 2 + ^CW, Pd m = \± 2 ~ ^CW- (13) 

The evolution equation for the energy density of the scalar field is obtained by taking the time 
derivative of its energy density (|12[h which, after some rearrangements, yields 

p + 3H(p <p +p c/> ) = - (p dm - Pdm) ^ - ' ( 14 ) 
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thanks to (|10|) and 1)12)1 ■ Following the same methodology, the time derivative of the energy density 
of the dark matter (|13(l i leads to the evolution equation for dark matter density: 

Pdm + 3H (p dm + p dm ) = (p dm - p dm ) i (15) 

9 

by using equations (|llfl and (|13|l 9. From equations l|14|) and l|15|) we infer that there exists a transfer 
of energy from the scalar field to the dark matter field which is given by the term (p dm — p dm )§ j '(f>. 
Moreover, the term 2^<fxpp/(l — £</> 2 ) in 1)14)1 is the responsible for the energy transfer from the 
scalar field to the gravitational field. 

One can obtain from equation 1)15)1 together with the definition p dm = m dm n dm = (,4> n dm a 
general expression for the particle number density for dark matter which reads 

n drn + 3Hn dm = -E^L I 2,H + - ] . (16) 
m dm y cj) J 

The above equation clearly simplifies to the usual n dm + 3Hn dm = for the vanishing dark matter 
pressure case and otherwise leads to a non conserved dark matter particle number density due to the 
particle creation by thermodynamic processes. In the next sections we shall find the cosmological 
solutions for this model by considering two cases, namely, a prcssurcless dark matter and a dark 
matter field with non- vanishing pressure. 

The conservation of the energy-momentum tensor of baryons and radiation with pb = and 
Pr = Pr/3 lead to the well-known relationships pb cx 1/a 3 and p r cx 1/a 4 . 



3 Pressureless dark matter 

For the case of a pressureless dark matter we have p dm = and it follows from equation (|16f) that 
p dm tx 4>/a? . Moreover, this last relationship implies that the particle number density of the dark 
matter is conserved, i.e., h dm + 3Hn dm = 0. Here we have only one equation to solve - namely, the 
evolution equation for the scalar field 1)10)1 - in order to obtain the time evolution of the acceleration 
field and of the energy densities of each component. By changing the variables, we can express the 
equation for the dark energy scalar field 1)10)1 in terms of the red-shift z, since a = 1/ (1 + z) by 
assuming that a(0) = 1. Hence it follows 

(1+^) 2 P ^ ^ + Z){P-ZP) ± , , dV , £0(p-3 P ) , pdrn =Q 



3(1 - ^0 2 ) 6(1- £02) # l-£0 2 

where hereafter the prime denotes the differentiation with respect to z. 

In order to solve ()17)l we have to choose a form for the potential density. Among the several 
existent models (see, for example |22]), we select a potential that was also used in the references |20l 
El El: 

where K is a constant. This form has the property of blowing up for small values of <f>, preventing 
the field - and therefore the dark matter energy density - from becoming negative. Furthermore, 
we have to specify initial conditions for the scalar field and its derivative as well as for the energy 
densities of the baryons, radiation and dark matter. The initial conditions for the energy densities 
were chosen from the present known values given in the literature (see j2U for a review) for the 
density parameters fi*(z) = Pl (z)/p(z), i.e., fl b (0) = 0.04995, Ct r (0) = 5 x 10" 5 and fl dm (0) = 0.23. 
The constant K was determined from the present value for the dark energy density parameter, i.e., 
O^(0) = 0.72, since it was considered that the present value of the scalar field has a small value, 
so that p<#)(0) ~ V. There still remain a and £ as free parameters, whose influence will be studied 
later in this section. 



4 



We have plotted in figure Q a phase-portrait that exhibits the evolution of a set of initial 
conditions. By varying the values of the scalar field in order to produce small variations in the 
initial dark energy density, it follows that some solutions evolve into a non-physical behavior since 
they lead to negative values of p^, whereas there exists another set of solutions that causes the dark 
energy density to increase too fast and is the responsible for almost all contributions to H and p. In 
order to get a physical acceptable behavior, it is necessary to impose a very small and positive value 
for the scalar field slope at early epochs, i.e., z>1000. Physically this means that the scalar field 
starts to roll more significantly around that time and then evolves until today. Moreover, negative 
derivative of the scalar field at high red-shifts would imply on negative dark energy density, due to 

03- 




Hubble parameter (H) 

Figure 1: Typical phase-portrait (p$, H) corresponding to the case a = 1/3 and £ = —0.2 showing 
the evolution for both physical (p^ > 0) and non-physical (p^ < 0) solutions. 

The red-shift evolution of the density parameters is plotted in the left frame of figure [21 We 
have considered a fixed value for the exponent a, namely a — 1/3, and different values for the 
coupling constant £. Two different cases were analyzed: (a) a slightly varying scalar field where the 
scalar field does not change significantly with the red-shift, showing that the differences from the 
usual quintessence model (£ = 0) are basically due to the influence of the non-minimal coupling 
(£ = —0.2); (b) a strongly varying scalar field, characterized by £ = —0.3, where the initial value 
of the scalar field has a smaller value than that of the former case but it has a more accentuated 
change when the red-shift increases. We infer from this figure that by increasing the red-shift the 
density parameter of the dark energy decays more slowly for the strongly varying case (£ = —0.3) 
followed by the slightly varying case when £ = —0.2 and £ = 0, whereas the density parameter of 
the dark matter increases more slowly for these cases. This can be understood by noting that the 
last term of the right-hand side of equation (|14(l - that depends on the coupling constant - reduces 
the energy transfer from the scalar to the dark matter field. For £ — the dark energy-dark matter 
equality occurs at z « 0.4, the dark energy-baryons equality at z f=a 1.4, whereas for £ = —0.2 
and £ = —0.3 these equalities occur at higher red-shifts. Moreover, for £ = the radiation-matter 
(baryons plus dark matter) equality occurs at z « 3700, whereas for £ = —0.2 and £ = —0.3 this 
equality happens for lower and higher red-shift values, respectively. 

The deceleration parameter q = —da /a 2 as a function of the red-shift z is plotted in the right 
frame of figure|21 We note that by increasing the coupling constant the transition from a decelerated 
to an accelerated phase occurs at lower red-shifts. According to present experimental values [5] 
the decelerated-accelerated transition occurs at zt = 0.46 ± 0.13 and we can observe from figure 
13 that for the slightly varying case with £ = —0.2 this transition is at zt = 0.46, whereas for the 
strongly varying case with £ = —0.3 at zt — 0.43. However, in both cases the present value of the 
deceleration parameter has a smaller modulus than the one given in the literature for the AC DM 
case, i.e., go ~ —0.55. 
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Figure 2: Left: evolution of density parameters vs. red-shift z, showing that the coupling 
changes dark energy-dark matter equality and radiation-matter equality with the red-shift; right: 
deceleration parameter q vs. red-shift z, the transition red-shift is closer to as the coupling 
increases. 



By changing the values of the potential exponent between the range < a < 5 it is not possible 
to see any significant difference in the behavior of the curves in figure |2 however, for larger values 
(a « 10), there is a slower decay of the dark energy density and a later transition to the decelerated 
phase. This occurs because V is of the same order as the other terms for red-shifts up to z « 4, 
and in this range the potential is mainly controlled by the value of K for a slightly varying <f>. 

In figureEOwe have plotted the difference /xo between the apparent magnitude m and the absolute 
magnitude M of a source, given by fj,o = m — M = 51ogc?i + 25 with 

d L = (l + z)cH^ / — - , (19) 
Jo H(z) 

as a function of the red-shift z, where d,L is the luminosity distance given in Mpc. In this figure 
the circles represent the experimental values taken from the work by Riess et al. [5] for 185 data 
points of super-novae of type la, whereas the dash-dot line stands for the standard KCDM model. 
We infer from the figure that all cases fit well at low red-shifts, but only the strongly varying case 
with £ = —0.3 stays near the KCDM over all the experimental range. 

To summarize, the best solutions are those that allow the scalar field to change only until the 
dark energy density parameter becomes small and the scalar field slope tends to a small positive 
quantity. 



4 Dark Matter with non-vanishing pressure 

Let us analyze the case where the dark matter has a non- vanishing pressure. For this end we recall 
that in the presence of matter creation the first law of thermodynamics for adiabatic (dQ = 0) open 
systems reads 

dQ = d{pV) + pdV - Pdm+Pdm d{n dm V) = 0, (20) 

where V denotes the volume and it was supposed that only dark matter creation is allowed. By 
considering V cx a 3 and the conservation laws for baryons and radiation it follows from (12011 

P,p + Pdm + 3i/(/50 + pdrn + P<t, + Pdm) - /3dm + Pdm (n dm + 3Hn dm ) = 0. (21) 

Hdrn 
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Figure 3: Difference fiQ vs. red-shift z, showing that coupling can increase or decrease the difference 
between magnitudes. 



If we compare the above equation with the sum of equations (|14|l and (|15(l , we obtain 



1 - £0 2 n d n 



'-{h dm + 'iHn dm ). (22) 



Now by using again (|15|l and p dm — m dm n dm — C4>n dm ^ we get the following expression for the 
pressure of the dark matter 




2 ^ (23) 



t<P)(3H + <j>/<t>) Pdm 

We call attention to the fact that only the plus sign must be chosen in the above equation in order 
to obtain a positive dark matter pressure. Moreover, we observe that in this case, when £ = 
we have from equation 12 I'll that the particle number density of the dark matter is conserved, i.e., 
"h-dm + 3Hn dm = and from (|23|l that the dark matter pressure vanishes. For the case analyzed in 
previous section there exists no constraint in the values of £ when we assume from the beginning 
that the dark matter pressure vanishes. 

We have now a system of coupled differential equations to solve, namely 

( 1 + 2 ) 2 P x» ( 1 + Z )(P-^P) x' , dV , ^(P-ZP) , Pdm ~ Pdm = Q , 24 v 



3(1 -^ 2 ) 6(1 -£0 2 ) r dcf> l-^ 2 



Hp dm H~ Pdm) ( \ y /or:\ 

Pdm 77- x = {Pdm ~ Pdm) "J, 25 

[1 + Z) (j> 

where p dm is given by 123|) with the positive sign. The initial conditions used here are the same as 
those of the previous section. 

Figure ^corresponds to figure[2]for the case of a non- vanishing dark matter pressure, where we 
have analyzed two cases (£ = ±0.05) and compared them with the minimally coupled pressureless 
case. We observe from the left frame of this figure that the decay of the density parameter of the 
scalar field for the negative coupling constant is slower than the one for the pressureless case, whereas 
for the positive coupling constant the decay is more accentuated. This can also be understood by 
analyzing the term — 2£(/></>p/(l — ^tfi 2 ) on the right-hand side of evolution equation for the energy 
density of the scalar field (|14fl . since it changes its sign when £ is positive or negative and it is 
equal to zero for the pressureless case. However, for both cases the radiation-matter (baryons plus 
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Figure 4: Left: evolution of density parameters vs. red-shift z; right: deceleration parameter q 
vs. red-shift z. The uncoupled case represents an intermediate case between positive and negative 
couplings for the density parameters at lower red-shifts (z < 2), and the coupling leads to a more 
recent transition redshift. 



dark matter) equality occurs for higher red-shifts than that of the pressureless case. From the 
right frame of figure 0] we infer that the transition from a decelerated to an accelerated phase for 
£ = —0.05 occurs at zt = 0.47 and the present deceleration parameter for this case approaches the 
value qo ss —0.55. 



5= -0.05 

5 = +0.05 




redshift (z) 



Figure 5: Difference vs. red-shift z where the coupling may increase the difference between 
magnitudes. 

In figure [5] it is plotted the difference between the apparent magnitude and the absolute magni- 
tude of a source as a function of the red-shift for the case of a non- vanishing dark matter pressure. 
We observe that the £ = +0.05 case stays near the ACDM model, whereas the £ = —0.05 case 
improves the behavior of the theoretical curve with respect to the experimental values for low 
red-shifts. 

Finally, the left frame of figure [B] shows the evolution of the dark matter pressure for the 
situations discussed above. We note that for £ = —0.05 the dark matter pressure vanishes when the 
red-shift goes to zero and tends to a small value for higher red-shifts. We emphasize that although 
low, the dark matter pressure introduced here has some significant influence on the behavior of the 
evolution of the whole system. For completeness we have plotted in the right frame of this figure 
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5=0 




redshift (z) redshift (z) 

Figure 6: Left: dark matter pressure Pdm/ Pdm vs. red-shift z\ right: dark matter particle number 
density ridm(z) I Vidm(0) vs. red-shift z, showing the influence of the coupling into the pressure and 
particle number density of the dark matter. 



the evolution of the dark matter particle number density as a function of the red-shift, and, as was 
expected, the particle number densities for the non-vanishing pressure cases increase more rapidly 
with the red-shift than that of the pressureless case. 

5 Final remarks and conclusions 

An equivalent description for the model studied in previous sections is to write Einstein field 
equations as 

Rpv — 2^9^ = ^liu ■ (26) 

In this case, the energy-momentum tensor of the sources is conserved, i.e., V^T^" = 0, whereas the 
energy-momentum tensor of the scalar field reads 



T+ = d^d^ - Qd CT <W - Vj g,, - £ (V M V^ 2 - g^VrVtf) + - \r 9 ^ ■ (27) 

By considering the Robertson- Walker metric and the representation for the energy-momentum 
tensor (T M „) = diag(p, — p, — p, — p) it follows the Friedmann and acceleration equations in their 
usual forms, i.e., 

H* = % a - = -e±&. (28) 

However, the energy density and the pressure of the scalar field in this description have new con- 
tributions that are given by the underlined terms: 

^ = ^|V + 6WI^, (29) 

P4> = \j> 2 - V - 2£ (# + 4> 2 + + jtfp , (30) 

so that the sums p = p^ + p drn + p b + p r and p = p^ + Pdm + Pb + Pr give 

- _ P<p + Pdm + Pb + Pr P ~ _ P4> + Pdm + Pb + Pr P 

showing the equivalence between this description and the former one. 
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To summarize, we first observe that this model has a strong dependence on the initial conditions. 
For the pressureless dark matter case we have shown that: (a) a stronger coupling between the scalar 
field and the scalar curvature leads to a decelerated-accelerated transition and matter-radiation 
equality at lower red-shifts; (b) a stronger scalar field variation implies that stronger couplings 
adjust better with the experimental data of luminosity distance, but the matter-radiation equality 
occurs at higher red-shifts. For the case of a non- vanishing dark matter pressure we have shown that 
negative small values of the coupling constant leads to: (a) a decelerated-accelerated transition at 
lower red-shifts with a better adjustment of the present value of the deceleration parameter, but the 
matter-radiation equality occurs at higher reds-shifts; (b) a better adjustment with the experimental 
data of luminosity distance at low red-shifts. 
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